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We revisit the problem of entanglement of thermal equilibrium states of composite systems. We 
introduce characteristic, viz. critical, temperatures -and bounds for them marking transitions from 
entanglement to separability. We present examples for the various possible thermal entanglement 
scenarios in bipartite qubit/qubit and qubit/qutrit systems. 
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I. INTRODUCTION AND SOME GENERAL RESULTS 

Consider a finite composite quantum system described by the complex Hilbert space 7i — Hi <8> H.2 8> ■ • • <£> H-n, 
which is the tensor product of N (> 2) finite dimensional Hilbert spaces Tij of dimension dj > 2 (J = 1, 2, • • • , N); 
and has dimension D = d\di ■ ■ ■ dw- Given any Hamiltonian H = H* acting on Ti, the thermal equilibrium (or Gibbs) 
state pr for temperature T is pr = exp(—H/T)/tr(exp(—H/T)). Since the quantum system is assumed to be finite, 
negative temperatures are possible, but they will be disregarded here for reasons of economy although the effects 
described below are also present for negative temperatures. 

Since the system has components, the question of the entanglement of the Gibbs-states arises. We recall that an 
arbitrary state uj (mixed or pure) of the composite system is said to be separable or unentangled if it can be written 
as a mixture (finite convex sum) of pure product states. If this is not the case, the state is said to be entangled. 
At present there are finite algorithms deciding whether a given state is entangled or not only for two qubits (N = 2 
with d\ — di =2); and for N = 2 with d\ = 2 and di — 3. The entanglement studies presented in the literature 
have, perforce, been restricted to investigations of the entanglement of the Gibbs-state reduced to the possible two- 
component subsystems, or to the study of part icular so called entanglement witnesses or other entanglement monotone 
functions [Jl3yiSHBSIaSll3lllll3ll3illJ|- f n section II, we present some features of the entanglement issue for 
thermal equilibrium states for the bipartite systems mentioned. The rest of this introduction collects basic information 
which is generally valid, i.e., for any N and D, and some of which seems to have been overlooked previously. 
In what follows we often identify the state to viewed as a linear positive functional acting on operators A and giving 
their expected- value u>(A), with the density operator u> such that ui{A) — tr(ujA). If the Hamiltonian is a real multiple 
of the identity, H = c-1, then px — 1/D for all T, or viewed as a functional, the normalized trace which we denote by 
r. Since 1/D = (1/di) £g> (l/cfe ) ® • ■ • ® the normalized trace r is a separable state. We assume henceforth that 

the Hamiltonian H is not a multiple of the identity, and let P_ be the spectral projection associated with the minimal 
(ground-state energy) eigenvalue s_(i?) whose multiplicity we denote by m_. For the ground-state po = P-/m-, one 
has limr_>o Pt = Po where the limit can be taken in various ways. As the limit of matrix elements in any orthonormal 
basis you wish; or limr^o t r (\ Pt — Po \) = where | X | denotes the absolute value of the operator X; etc. Moreover 
limr^oo Pt = Poo = t in the same sense, and the map < T i— ► pr is continuous. 

For what follows it will be important to observe that the energy U(T) = Pt{H) as a function of temperature, is a 
monotone increasing continuous function with limT^o U(T) — U(0) — S-(H) and limx^oc U(T) = t(H) = tr(H)/D. 

A. Various "critical" temperatures 

A critical or threshold temperature above which entanglement is impossible has been observed in all cited studies. It 
is shown in 0] that for every Hamiltonian there exists a finite critical temperature which satisfies: (i) < T5 < 00; 
(ii) pr is separable for every T € [T$, 00]; and (iii) for every < T < Ts there are entangled thermal states with 
temperatures in the interval [T,Tg). We call Ts the lower separability temperature. Even if H is not a multiple of 
the identity, it could still be trivial in the sense that it contains no interactions whatsoever between the component 
subsystems, H — J2n=i H^ n \ acting only on the n-th component. In this case px = Pj^ ® Pt^ ® * " ® Pt^ ^ s a 
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product-state, hence unentangled and Ts = 0. A very poor upper bound on Ts was given in |15j : and the question of 
what happens below Ts was posed and left open. It will be essentially answered in section II. 

It is important to stress the fact that the convex set of separable states is closed, and this implies that the set of 
non-negative temperatures T for which pt is unentangled, is closed. Accordingly, the non-negative temperatures T 
for which px is entangled is an open set in [0, Ts). 

By their very definition, entanglement witnesses or entanglement monotone functions will always provide temperatures 
for which the thermal states are entangled and these temperatures are lower bounds on Tg. A particularly simple 
witness is the energy itself as observed in Let rj = ini{u>(H) : lo is separable}, that is the lowest energy 

expectation value obtainable with an unentangled state. The infimum is assumed, and since the map to 1— > w(H) is 
convex-linear, it can be taken over the pure product-states. 

If for some T\ > we have pt x (H) < 77, then p Tl is entangled. By the monotone increase of U(T) and the intermediate- 
value theorem there is a unique Th > Ti such that U{Th) — rj and all Gibbs states with temperatures in [0, Th) are 
entangled (in particular p is entangled). In E2, Th is denoted by Te and called the entanglement-gap temperature. It 
follows that Th < Ts. It does not follow that Gibbs states with temperatures (immediately) above Th are necessarily 
separable. Although Th > does indeed signal the presence of thermal entanglement at low enough temperatures, 
its importance should not be overrated. The correct critical value is: Te = inf{T > : pt is separable}; for which 
one can prove, as in EH, that: (i) < Te < Tg, and pt e is separable; and (ii) Te > if and only if p is entangled, 
and in this case all thermal states with temperatures in [0, Te) are entangled. Alternatively, Te could be defined as 
the greatest temperature such that all Gibbs-states with temperatures below it are entangled; we call Te the upper 
entanglement temperature. Obviously, Th < Te, but one should expect that, in general, Th can be a rather poor 
lower bound on Te. The following example should serve as illustration. 

Consider two qubits; and suppose the minimal energy s_ (H ) of your Hamiltonian is doubly degenerate with ground- 
state vectors ipi = a® a, and 4>2 = (a ® + P ® ct)/V% where a (resp. (3) is an eigenstate of 173 to the eigenvalue 1 
(resp. —1) for one qubit. Then the ground-state is p = (1/2) | ipi){^\ | +(1/2) | ip 2 ) (^2 I, an d it is entangled (the 
partial transpose of po has a negative eigenvalue); rj = S-(H) — {if!i,Hipi) and thus Th = 0; but Te > 0. 

B. The "more mixed than" ordering of thermal states. Upper bounds on Ts 

The theory of the "more mixed than" partial ordering of states of a quantum system, [l6j . can be put to use 
in the discussion of entanglement of Gibbs states. It is a result of Wehrl and Uhlmann (cf. Refs. [l6l llT\). that 
< T < T' < 00 implies F(px) < F (pt'), for every unitarily invariant, concave, continuous real- valued functional F 
defined on states. It is shown in [l5j, that for every such functional F, for which F(us) — F(t) implies lo = r, there 
is a critical constant Cf < F(t) such that: (i) If the state to satisfies F(u>) > Cf then u> is separable; and (ii) For 
every possible value C of F below Cf there is an entangled state <j> with F((f>) = C. There is an analogous version 
of this for unitarily invariant, convex, continuous real-valued functionals. Thus, every unitarily invariant, continuous 
real-valued functional which isolates t and is either convex or concave, acts as a separability detector and can be 
used to obtain an upper bound on Ts. Indeed, take such a concave separability detector F. Then, T 1— > F(pr) is a 
non-decreasing continuous function for which limT^oo F(px) — F(t) > Cf- By the intermediate- value theorem there 
is Tf < 00 such that F(pt f ) — Cf and all Gibbs states with temperatures in [Tp,oo] are unentangled. It follows 
that T s <T F . 



C. Thermal entanglement scenarios 

Many of the available studies, e.g. 0, H 0, IE EH El, observe that the bipartite subsystem entanglement of 
thermal states need not be monotone in temperature. 

The above results allow one to distinguish various entanglement scenarios. The uninteresting scenario occurs when 
Te = Ts = as happens when the Hamiltonians present no interactions whatsoever. But this scenario is possible 
even when interactions are present, cf. Section II. 

The next scenario is that in which the ground-state is separable, i.e., Te — 0, but Ts > 0. Then as temperature 
increases away from zero, separability is lost at some temperature < T± < Ts. For temperatures in the "separable 
segment" [0,Ti] all Gibbs states are unentangled, and the segment (Ti,Ts) contains temperatures for which 
the corresponding thermal states are entangled. The possibility arises for various closed "separable segments" 
([0,Ti], [T 2 ,T 3 ], ■ • • [T n , T n+ i]) alternating with open "entanglement segments" ((Ti,T 2 ), ■ ■ ■ (T n+l ,T s )). We will 
present examples for this "abnormal" scenario in section II. 



3 



The other scenarios occur when the ground-state po is entangled, i.e., Te > 0. The normal case is Te = Tg and 
this is what has been observed in most studies we know of. Although, entanglement usually decreases with increasing 
temperature, one can have non- monotonous behaviour, cf. section II. The abnormal scenario is < Te < T$, 
and then there is a temperature T\ with Te < T\ < T$ such that for T £ [Tg,Ti], p T is unentangled but there 
are temperatures T" <E (Ti,Tg) for which px> is entangled. Again the way is open for closed separable segments 
alternating with open entanglement segments. Examples of this behaviour are given in section II. 



D. The modulus of separability 

Consider any state ui of the composite system and consider the segment joining the normalized trace r to lu, i.e., 
uj t = t ■ u> + (1 — t) ■ t with < t < 1. Since luq = r is separable, one will ask in case u>\ — lu is entangled, when 
as t increases does one lose separability? This question has been analyzed by many authors, notably by Zyczkowski, 
Horodecki, Sanpera, and Lewenstein. Il8l . who develop it to obtain a method of estimating the "size" of the separable 
states; and by Vidal , and Tarrach, who give a virtually complete treatment of the problem. In Ref. the 
modulus of separability of ui was defined as £(lu) = sup{0 < t < 1 : uif is separable }, whereas the quantity considered 
by Vidal and Tarrach is R(lu || r) = — 1 and called by them the random robustness of entanglement. Here 

we only need to observe that < £(lu) < 1, with £{lu) = 1 if and only if li is unentangled. Moreover, the upper- 
semicontinuity property of £ obtained in |15|, guarantees that the map T i— > £(pr) is continuous. This in turn, proves 
the claims made in £11 Al about the sets of temperatures where the Gibbs state is separable, respectively entangled. 

II. THERMAL ENTANGLEMENT IN QUBIT /QUBIT AND QUBIT/QUTRIT SYSTEMS 

Vidal and Tarrach, have computed the modulus of separability for a qubit/qubit (D = 4) or a qubit/qutrit 
(D = 6); their beautiful formula is 



v ' 1 + D | min{A(w),0} | ' 

where X(lu) is the minimal eigenvalue of the partial-transpose of lu with respect to the qubit. The plots which we will 
exhibit show T i— » £{pr) for selected Hamiltonians which exemplify the distinct scenarios. 

The critical values (or good upper and lower bounds) are known for some functional (e.g. partial ordered eigenvalue 
sums, von Neumann entropy) from Ref. and can thus be used to provide upper bounds on T5 as described in 
i llBI The map u> 1— > tr(uj 2 ) is strictly convex (cf. Ref. [l^. p. 47), unitarily invariant and is easily seen to isolate r. 
Furthermore, the critical value for the trace of the square is known to be 1/3 for two qubits and it lies between 1/5 
and 7/32 for a qubit/qutrit system (^H)- It turns out to be a rather useful separability detector because it is easy to 
calculate. Although it often provides a rather poor upper bound on T5, in all the examples to be presented below it 
gave better bounds than those obtained using the partial eigenvalue sums, or the von Neumann entropy. We denote 
by , the numerically obtained value of Tp -recall til Bh for F equals minus the trace of the square (using the bound 
1/5 for the qubit/qutrit case). We have not invested a lot of effort in the computation of 77 of £11 Al The temperature 
value Th obtained from the numerical, possibly inexact, value of 77 need not necessarily be a lower bound for Te or 
even Tg; but in our experiments this never happened. The eigenvalues of the Hamiltonian counting multiplicities are 
given as a row- vector h. Since in our definition of px we have incorporated Boltzmann's constant in the Hamiltonian, 
the components of h have the same dimension as the temperature. Since thermal equilibrium states are invariant 
with respect to addition of a multiple of the identity to the Hamiltonian, we choose S-(H) = 0. The eigenvector to 
the j-th eigenvalue hj is listed as a row vector ej , where the coordinates are with respect to the canonical orthonormal 
tensor-product basis built from the orthonormal basis {(1,0), (0, 1)} of C 2 , and {(1,0,0), (0, 1,0), (0,0, 1)} of C 3 . 

The Hamiltonians to be presented are specifically chosen to exhibit transitions from entanglement to separability 
below Ts, that is in the interval [0,Ts). The general idea is, obviously, to choose the eigenvector associated with the 
non-degenerate ground-state energy to be either separable or entangled and then the eigenvector associated to the 
first excited state to be, correspondingly, either entangled or separable. 
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A . qubit / qubit 

The minimal value of the separability modulus for two qubits is 1/3. 
Figure 1, shows the = Te < Tg scenario with a separable segment [0,Ti] followed by an entanglement segment 
(Ti,Tg). Using the very same eigenvectors as those of Fig. 1, but with h = (0, 1.5, 2, 3), px is unentangled for every 
T > 0, that is Tg = 0. 

Figure 2 shows two instances of the standard normal scenario (T# <) Te = Tg. 

In the qubit/qubit system, we have not found the scenario where < Te < Tg (cf. Figure 3 for a qubit/qutrit). 

B. qubit/qutrit 

The minimal value of the separability modulus for a qubit/qutrit systems is 1/4. 
Figure 3 shows the < Te < Tg scenario; the entanglement segment [0, Te) is followed by a separability segment 
[Te, 2~i] and a second entanglement segment (Ti, Tg). In this an further examples not shown here, we have found that 
once a region of the h-space is found where the pertinent scenario is present, the values Te,Ti,Ts are quite stable 
with respect to changes in h which are large with respect the temperature values. The inset of Figure 3 shows this 
effect. Moreover, the same eigenvectors used in Fig. 3, but with h = (1.7, 0, 1.75, 2, 3, 4) give Te — Tg — 0.699. 
Figure 4 shows the = Te < Tg scenario (cf. Figure 1 for the two qubit system) with four transitions: [0,Ti], and 
[T2,T 3 ] are separable segments alternating with entanglement segments (T l7 T 2 ) and (T 3 ,T S ). 



III. CONCLUDING REMARKS 

We have introduced two characteristic temperatures Te and Tg which organize the entanglement behaviour of the 
thermal state associated to any Hamiltonian of an arbitrary composite system. For qubit/qubit and qubit/qutrit 
systems, we have exemplified the possibility of various transitions from entanglement to separability as temperature 
increases from zero to Tg, above which entanglement is impossible. One could expect that the features found here 
will persist and be enhanced as N or D increase (although at present there is no manageable criterion to decide when 
a given state is entangled or not). Thus, in general, there will be many separable temperature segments alternating 
with entangled ones for multipartite systems of higher dimensions. 

The main motivations for the present study came from certain problems posed in Ref. |l5j. From that point 
of view, the main conclusion to be drawn from our findings here are theoretical and concern the results briefly 
mentioned in Ml Bl It was asked in Ref. |l5|: given a separable state lo which is not pure, does there exist a unitarily 
invariant, concave, continuous real-valued functional F defined for the states of the composite system which isolates 
the trace and such that F(lo) > Cf 1 - The answer given here is definitely no! Take any unentangled thermal state 
PTi such that for some T 2 > T\ the Gibbs state pt 2 is en tang led. Then there cannot exist an F with F(pr 1 ) > Cf 
because by Wehrl's result (cf. Ml F(pt 2 ) > Cf and by the separability of pt 2 would follow. 
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FIG. 2: Solid line: ei = (1/^2,0,0,1/^2), e 2 = (1/^,0,0,-1/^2), e 3 = (0,1/^2,1/^,0), e 4 = (0,1/^2,-1/^2,0); 
h = (0.75,0,0.75,2). T H = 0.73, T E = T s = 0.97, and T, = 1.04. Dashed line: ei = (1,0,0,0), e 2 = (0, l/y/2, l/y/2, 0), 
e 3 = (0, 1/72, -1/V2, 0), e 4 = (0, 0, 0, 1); h = (0.01, 2, 0, 4). T H = 0.377, T E = T s = 1.823, and T* = 2.181 
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FIG. 3: ei = (1,0,0,0,0,0), e 2 = (0, 0, x, 0, x, y), e 3 = (0,0,1/^,0,-1^/2,0), e 4 = (0 ,1/72,0,1/^2,0,0), e 5 = 
(0,l/\/2,0,-l/V^,0,0), e 6 = (0,0,j//V2,0,j//V2,-a;V2), where x = 0.2 and y = s/l - 2x 2 ; h = (0.75,0,0.75,2,3,4). 
Th = 0.13, Te — 0.296, Ti — 0.334, Ts = 0.571, and T« = 2.76. Inset: using the same eigenvectors and eigenvalues 
except, from top to bottom, hi = 0.75, 1 and 1.5 
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FIG. 4: ei = (1,0,0,0,0,0), e 2 = §(0,1,0,1,1,1), e 3 
e 6 = i(0, -1,0, 1,1,-1); h = (0,0.7,7,0.9,1,1.5). T H 
T, = 2.645. 



= (0,0,1,0,0,0), e 4 = i(0, 1,0, 1,-1,-1), e 6 = ±(0, 1, 0, -1, 1, -1), 
= T E = 0, Ti = 0.0355, T 2 = 0.467, T 3 = 0.476, T s = 0.923, and 



